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Abstract We compare three methods to calculate the nucleon-nucleon t- 
matrix based on the three-dimensional formulation of Ref. [1]. In the first 
place we solve a system of complex linear inhomogeneous equations directly 
for the t-matrix. Our second method is based on iterations and a vari- 
ant of the Lanczos algorithm. In the third case we obtain the t-matrix in 
two steps, solving a system of real linear equations for the k-matrix expan- 
sion coefficients and then solving an on-shell equation, which connects the 
scalar coefficients of the k- and t-matrices. A very good agreement among 
the three methods is demonstrated for selected nucleon-nucleon scattering 
observables using a chiral next-to-next-to-leading-order neutron-proton po- 
tential. We also apply our three-dimensional framework to the demanding 
problem of proton-proton scattering, using a corresponding version of the 
nucleon-nucleon potential and supplementing it with the (screened) Coulomb 
force, taken also in the three-dimensional form. We show converged results 
for two different screening functions and find a very good agreement with 
C*~) . other methods dealing with proton-proton scattering. 
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1 Introduction 

The nucleon- nucleoli (NN) t-matrix appears not only in the description of the 
NN scattering process but constitutes a key ingredient of the three- and four- 
nucleon calculations (see for example [HIS]). No wonder that many methods to 
calculate this quantity have been proposed and realized both in coordinate 
and in momentum space. For many years it has been natural to use the 
rotational invariance of the NN force and introduce a partial wave basis. 
This procedure has a clear physical meaning and is very accurate at low 
energies. However, at higher energies very many partial waves are necessary 
to achieve converged results, so methods employing vector variables might 
become more efficient. 

There have been several approaches to NN scattering avoiding partial 
wave decomposition (PWD). Pioneering calculations of Ref. [J] carried out 
for different Malfliet-Tjon-type potentials gave interesting numerical insight 
in the properties of the three-dimensional (3D) t-matrix and paved the way 
for further investigations including spin degrees of freedom. A helicity for- 
malism directly linked to the total NN spin was proposed in [5]. Later, the 
spectator equation for relativistic NN scattering was solved in [5] also using 
a helicity formulation. Closely related are 3D formulations for pion-nucleon 
scattering [7] and for scattering of protons on light nuclei [5] . 

These developments were then used in few-body calculations. The Fad- 
deev equation for the system of three bosons interacting via scalar forces was 
solved for the bound state case [5] as well as for three-body scattering [TU] 
directly in three dimensions. The helicity formalism of Ref. [5] was employed 
in 3N calculations in [TTirT^l[T5] . Recently we proposed a 3D framework for 
3N bound [14UT5] and 3N scattering states [16] , where spin-momentum oper- 
ators are introduced and treated analytically so the Faddeev equations turn 
into a finite set of coupled equations for scalar functions depending only on 
momentum vectors. For this formulation it is crucial that the most general 
form of the NN interaction can only depend on six linearly independent spin- 
momentum operators. In [I] we formulated a new approach to NN scattering 
and provided its numerical realization based on the matrix inversion method. 
A comparison with the standard PWD calculations for two quite different NN 
potentials and a few projectile energies proved that the new scheme is very 
accurate pQ. 

However, in view of expected applications to the 3N calculations, where 
we need the full off-shell t-matrix, we decided to develop a more efficient, 
iteration approach to calculate the 3D t-matrix. We also realized an old idea 
presented, for example, in [17; and prepared a scheme for obtaining the on- 
shell t-matrix from the solution of the k-matrix equation. 

In Sect. [5] we repeat briefly the main points of our approach starting from 
the most general form of the NN potential. Next we derive the correspond- 
ing equation for the k-matrix and formulate the subsequent equation for the 
on-shell t-matrix. Finally, we formulate our iteration scheme. Numerical real- 
izations of our different approaches that employ a recent chiral next-to-next- 
leading order (NNLO) NN force [T81I19U20) are presented in Sect. H Here, 
in addition, we apply one of our schemes to the problem of proton-proton 
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scattering (including the strong and screened Coulomb forces), and compare 
results for two different screening functions. We conclude in Sect. [U 



2 The Formal Structure 

For the convenience of the reader we repeat here the main points of our 
formalism, introduced in Ref. [T]. The NN potential in the two-nucleon (2N) 
isospin space, spanned by the four states | tm t ) (| 00), | 1 — 1), | 10), | 11)), 
is given as 

{t'm t , \ V\tm t } = 6 tt >6 mtmtl V tmt , (1) 

allowing for charge independence and charge symmetry breaking. The most 
general rotational, parity and time reversal invariant isospin projected NN 
force is then expanded into six scalar spin- momentum operators w;(<xi, er 2 , p' , p) 
[2] as 

6 



V^=5j«r t (P'.P)«'j(ffi.o-a,p / ,p), (2) 



using [1 



■w 1 {a 1 ,<T2,p',p) = 1 , 

l0 2 (o-l,<T 2 ,p',p) = (Ti ■ <T 2 , 

?« 3 (cr 1 ,(T2,p / ,p) = i (<7i +er 2 ) ■ (p x p') , 
W4(eri,er 2 ,p',p) = rr\ ■ (p x p') cr 2 • (p x p') , 
Ws(0-1,0-2,P',P) = CTl • (P' + P) <*2 ' (P' + P) , 
w 6 (eri,cr 2 ,p',p) = <T\ ■ (p' - p) <7 2 • (p' - p) , (3) 

and scalar functions vp (p', p) which depend only on the p and p' momenta. 
A corresponding expansion into the iUj(<Xi, cr 2 , p', p) operators can be applied 
to the NN t-operator at the 2N internal energy z: 

6 

t tm t {z) = t *p ( P ', P ; z) Wj <r 2 , p', p) . (4) 

The z parameter appears in the Lippmann-Schwinger (LS) equation 

t tmt (z) = V tmt +V tmt G Q (z)t tmt (z), (5) 

with Gq(z) = (z — Ho)^ 1 being the free resolvent. 

Inserting Eqs. © and (HJ) into the LS equation ([S]), operating with Wk(fx, cr 2 , p', p) 
from the left and performing the trace in the NN spin space leads to 



E 



A fej (p' )P )tf V,p;z) =E^(P / 'PH* mt (P / .P) 



/ 



d 3 p" V T* (P'> P") tip (P". P; B kjf (P\ P", P) • (6) 
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The scalar coefficients Akj and Bkjy , denned as 
Bkjji (p',p",p) 

are given explicitly in Appendix A of Ref. [T]. In this way we obtain a set 
of six coupled equations for the scalar functions i^ mt (p', p; z), which depend, 
for fixed z and |p|, on two variables, |p'| and the cosine of the relative angle 
between the vectors p' and p, given by p' • p. 

In Ref. [T] we showed all steps leading to the NN elastic scattering ob- 
servables. The LS equation was solved there for positive energy of the NN 

2 

system, z = E c . m . = where m is the nucleon mass. We chose the following 
representation for the vectors p, p' and p": 

p= (0,0,1), 
p' = (Vl-x' 2 ,Q,x'), 

p" = (y/l-x" 2 cosip", y/l-x" 2 smtp",x"), (9) 

and prepared grids for the different variables in the problem. We took n x 
Gaussian points for the x" integration and used the same grid for the x 1 
points. We employed n p Gaussian points for the p' and p" grids, which are de- 
fined in the interval (0,p). In addition, po is added to the set of p' points. The 
intermediate ip" integration was performed with n v >> Gaussian points. Thus, 
from Eq. ([B]) we got a system of 6 x (n p + 1) X n x linear equations for given 
Po and fixed p. In Ref. [1] we set from the very beginning p = po; so the solu- 
tion contained the on-shell t-matrix in the operator form, tj(po,po, x'; E c , m ). 
Since the six operators © are on-shell linearly dependent, our solutions were 
not unique. They always led, however, to stable and unique predictions for 
the observables. In Ref. Q] we worked with the standard LU decomposition of 
Numerical Recipes [22] for two different NN forces (the Bonn B [23] potential 
and a chiral NNLO potential [18, 19, 20 J and made calculations for several 
E c . m . energies. For all considered cases we obtained a perfect agreement be- 
tween this new 3D approach and the calculations based on standard partial 
wave methods. 

Solving Eq. ([5]) directly as a system of inhomogeneous coupled algebraic 
equations is not efficient. Especially in view of the applications in the three- 
nucleon system, where a full off-shell t-matrix is needed, it is worth looking 
for a faster method. A typical job employing n x = 36, n p = 36 and n v " = 60 
Gaussian points might take a few hours on a PC. In Ref. [T] we tried to solve 
Eq. ([6]) by iterations. In the matrix form this equation was written as 



Tr(wk(<Ti,(T2,p',P) w J (cTi,cr2,p',pyj , (7) 
Tr (w k (a- 1, <r 2 , p', p) Wj((Tx,<T2, p', p") 

XW i /(<7i,tT 2 ,p",p)) , (8) 



At = Av + Bt , 



(10) 
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where tj(p' ,p,x' ; z) (Vj(p',p, x')) constituted the components of the t (v) 
vector and 

P 

{Bt) k {p' ,p,x')= / dp"p" 2 — -2 fk(p";p',p,x'), (11) 

J - p" + le 

with 

fk(p";p',P,x') = 

g 1 27T 

3J'=1_1 o 

(12) 

where 

y = x'x" + \fl - i'Vl - x" 2 cos ip" . (13) 

It was possible to use such p, p' and x' points that the matrix A could be 
inverted [1]. In this case Eq. (O was written as 

t(p',p,x')=v( P ',p,x')+A- 1 (p>,p,x') (Bt)(p',p,x') (14) 

so we arrived at the following iteration scheme: 

t^(p',p,x')=v{p',p,x'), 
tW(p',p,x')=v(p',p,x') 

+ A- 1 (pf,p,x , )(Bt( n - 1 ))(p',p,x'), forn>l. (15) 

We found it, however, very difficult to maintain numerical stability for this 
iterative method. Fortunately, it turned out that our scheme could be easily 
modified to give stable results. It is sufficient to combine, on the level of 
analytical calculations, the A k j and B k jj> coefficients, preparing 

6 

B kjj , (p', p", p) = *ia(p': P) B W (p'» p"> P) ■ ( 16 ) 
i=i 

The new, numerically safe, iteration scheme reads thus: 
tW(p',p,x')=v(p',p,x') 

t^(p',p,x') = v(p',p,x') + (M n -V)(p',p,x'), for n > 1, (17) 

where 

v 

(Bt) k (p',p,x>) = f dp"p" 2 - f k (p";p' )P ,x') (18) 

J Pq - p " + ie 
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and 

f k (p";p',p,x') = 

6 1 27r 

m J dx " J V B kjj , (p',p",p, x', x",cp") Vjtfy, y) t r (p",p, x") . 
j'.j'=i_i o 

(19) 

The above given iterative method proves also very fast. It is very important 
from the numerical point of view that the summation over j and the integral 
over if" do not involve tj(p',p,x'), which means that they can be prepared 
and stored in advance and later used in each iteration. Furthermore, this 
scheme can be easily implemented on a parallel machine and used to calculate 
tj(p',p, x') for many p points at one shot. In the case, when one is interested 
in the on-shell t-matrix, it is impossible to set directly p = p$, if the p' points 
contain pq. It is sufficient to make calculations for two p values, p — po ± S Po 
(S Po w 0.01 fm _1 ) and then take the average: 



tj(po,Po,x';z= — ) w i (t 3 (pa,p - 5 Po ,x';z = — ) 



1- 1 ' 

m ' 

+t j (p ,po + 5 Po ,x';z=^)j . (20) 

This method was successfully employed in the recent 3N bound state cal- 
culations [15] , In this case the kernel contains no pole, which makes the 
calculations even faster. We use a variant of the Lanczos method [H] to sum 
the resulting Neumann series, which proves to be more accurate than the 
standard Pade scheme. 

Alternatively, in many calculations using PWD of the NN potential, the 
on-shell t-matrix is obtained not directly but in two steps. First one defines 
the k-matrix through the principal value kernel |17] 

(p'|fc tm '(z)|p) = (p'|^ tmt |P> 

d 3 p" (p' | v tmt | p") — <P" I fc tmt (*) I P) (21) 

' m 

and the connection to the t-matrix is given as 
(p'\t tm '(z)\p) = (p'\k tm Hz)\p) 

-m J d 3 p" <p' | (z) | p") 6^z - ^ J (p" | t tm * (z) | p) . (22) 

2 

If we take z = — , |p"| = |p'| = |p| = p$, then we get the following (on shell) 
relation 

(p'\t tmt {z)\ V ) = (p'\k tmt {z)\ V ) 

- l -nmp Q J dp" <p' | fc tm ' (z) | p") (p" | t tm ' (z) | p) , (23) 
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where dp" denotes the two-dimensional angular integration. 

Including spin degrees of freedom and using the expansion of the k-matrix 
into the iu^o - !, <r 2 , p', p) operators (jSJ), it is easy to find the equation corre- 
sponding to (|2"TT) but for the scalar coefficients of the k-matrix, fc* m *(p', p; z): 

X> i (^ P )*f i V ) p;*) = E^y(p''P>rv>p) 

3 3 

+ J d 3 p" £ «f * (p', p") — ^ kf- (p", p; z) S W (p', p", p) . (24) 

jj' Z m 

The coefficients A^j and B^y are the same as in Eq. ([B]). By simple in- 
terpolation or by taking the average over two p values, we get the on-shell 

k-matrix coefficients, kj mt (po,po,x';^). In order to formulate an equation 
corresponding to (|2"5)) . we have to take into account that on-shell there are 
only five independent «?j(<xi, <T2, p', p) operators. This follows directly from 
the (on shell) relation 

CTi • er 2 = -, — ^ o-i • (p x p') a 2 ■ (p x p') 

IP x P ) 

+ 7 i M2 °"1 • (P + P') °"2 ' (P + P') 

(P + P') 

+ (p_ p /)2 °"i • (P - P') o-2 • (p - p') , (25) 

and can be written as 

1 1 1 

w 2 = — T7 + ^— tt; 77^5 + —^7- — U>6 . (26) 

Pt(l-x> 2 ) 2p 2 (l + x') 2p 2 (l-x') 

Using Eq. (1261) it is easy to write the on-shell k-matrix as a linear combination 
of only five operators. Then the on shell equation corresponding to (|23l) takes 
the form 

5 2 

Pi- 
rn 

3=1 



T / A kj(Po,Po,x')tf lt (Po,Po,x';^) = 

3=1 

5 2 
3=1 

1 2tt 5 



3>3'=1 

x ^(po.Po^'';^)^'^,^,^^',^^''). (27) 

In principle, any operator appearing in (1261) can be eliminated but not 
all the possible choices are numerically equivalent due to different behavior 
of the resulting coefficients in the vicinity of x' — ±1. This is demonstrated 
in Sect. [3l In the case when the operator or w§ is removed, the A^j and 
B).jji coefficients in Eq. ((27)) have to be renumbered accordingly. 
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3 Numerical Results 

In this section we focus on the on-shell t-matrix obtained without employ- 
ing PWD. We do not compare our 3D results with the results obtained by 
standard momentum space PWD, since this comparison was carried out suc- 
cessfully in pp. 

We show first that the methods discussed in the previous section lead 
to stable numerical results. We use the matrix inversion scheme, introduced 
already in [T] . In present calculations we use a parallel version of the matrix 
inversion, based on the SCALAPACK library [26]. We implement also the 
iterative algorithm and sum up the resulting Neumann series using a variant 
of the Lanczos method [23] . Finally, we solve the equation for the k- matrix 
and, going on-shell, eliminate either the wq or the W4 2N operator. 

We perform calculations using all these methods for two projectile lab- 
oratory kinetic energies, 13 and 300 MeV. We use first the neutron-proton 
version of a chiral NNLO NN potential [TgllTW!?U] . described briefly in Ap- 
pendix A of Ref. pp. 

In Figs. Q] and [5] we show results for several NN scattering observables. 
In the case of 13 MeV, the agreement among the four types of results is 
very good both for the neutron-neutron (left panel) and the neutron-proton 
(right panel) system. For the energy of 300 MeV we see that the calculations 
using the k-matrix equation without the wq operator on-shell, yield slightly 
different results for the R and A observables in the case of the neutron-proton 
system, than the other three schemes. These numerical instabilities become 
better visible for the Wolfenstein parameters |17] c, g, h and m shown in 
Fig. |21 Thus this particular numerical scheme should be avoided. 

In the case, where iterations are used, we wanted to know how many it- 
erations are needed to achieve full convergence of our results. We expected 
that the convergence depends on the NN system internal energy and show 
in Figs [J] and [5] our results for various numbers of iterations again for the 
small (13 MeV) and large (300 MeV) laboratory energy. At 13 MeV and for 
the neutron- neutron system we find that 10 iterations are fully sufficient. 
For some observables results obtained with 8 iterations are very close to the 
predictions based on 20 iterations. For the neutron-proton case one needs 
generally 12 iterations. At 300 MeV we see even faster convergence. For the 
neutron-neutron case predictions based on 6 iterations have already a correct 
shape (this was generally not the case for 13 MeV) and 8 iterations bring us 
very close to the fully converged results, although some tiny deviations are 
still visible. For the neutron-proton system results obtained from 10 itera- 
tions are still slightly different from the fully converged predictions. We can 
thus state that for energies below the pion production threshold, one needs 
no more than 12 iterations to arrive at fully converged results. Since the 
essential part of the iteration kernel can be prepared in advance, iterations 
run very fast, especially for the case of negative 2N internal energy, where 
no singularity is present. This method is considered by us as the best for 3N 
calculations. 

In the remaining part of this section we consider the most difficult case in 
the nucleon-nucleon scattering, namely the proton-proton case with inclusion 
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of the Coulomb force. As described in Ref. [27] , contrary to the approaches 
using the so-called renormalization procedure (see for example 0), we add a 
screened Coulomb potential to the (strong) proton-proton 2N force given in 
the operator form and use this sum as input in our 3D calculations. We work 
with the proton-proton version of the same chiral NNLO 2N potential and 
choose a relatively small laboratory energy (13 MeV), where Coulomb effects 
are expected to play an important role. We assent to two types of screening, 
the exponential screening, si(r;n, R), 

Sl (r;n,i?)=exp(-(^)") , (28) 

restricting ourselves to the case n = 4 and the localized (transition from 1 
to takes place in a finite interval) screening, S2(r; R), proposed in [5]: 

s 2 (r; R) = 0(R - r) + ±6(r - R)0(3R - r) (sin (|^) + l) , (29) 

with 6{x) being the (unit) step function. The 3D momentum space screened 
Coulomb potential is obtained as in Ref. [2"S] . 

In Figs. [5] and [7] we demonstrate a convergence of the results for some 
observables with the increasing parameter R. In the case of the exponential 
screening with n = 4 it is achieved for R= 120 fm, while for the second type 
of screening R= 55 fm is sufficient. In Fig. [5] we show that, as expected, 
the two screening prescriptions lead to the same final result for the selected 
proton-proton observables in our fully 3D calculations. The inset in Fig. [5] 
serves to compare the two screening functions for the parameters used in the 
actual calculations. 

Proton-proton scattering was already considered by us in j27) . There we 
used the Vincent-Phatak method [H] and the method combining PWD cal- 
culations with the full 3D t-matrix of the screened Coulomb potential. Using 
those two methods, we performed calculations for the proton-proton strong 
potential and screening functions used in this paper. Predictions for the R 
parameters yielding converged results in the case of the fully 3D calculations 
(solid line) are compared in FigsJS] (exponential screening) and ITU1 (localized 
screening) with the Vincent-Phatak method (dotted line) and the method 
combining PWD and 3D approaches (dashed line). We observe a good or 
even a very good agreement among all the methods for the polarization ob- 
servables, especially for the exponential screening. Discrepancies in all the 
cases, including the most visible ones for the cross section, do not exceed 1- 
2 %. This agreement provides yet another proof that our numerical treatment 
of nucleon-nucleon scattering is fully reliable. 

4 Summary and conclusions 

In pQ we formulated an approach to the NN system employing spin-momentum 
operators multiplied by scalar functions, which only depend on the momen- 
tum vectors. This representation of the NN potential leads to a system of 
six coupled equations for the scalar functions defining the NN t-matrix. This 
system of equations was first solved using the matrix inversion method. 
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In the present paper we formulate two other methods to calculate the 
nucleon-nucleon t- matrix in the same 3D formulation of Ref. pQ. Our sec- 
ond method is iterative and uses a variant of the Lanczos algorithm 24J. 
In addition, as a third possibility, we treat the equation for the k-matrix 
and obtain the on-shell t- matrix coefficients, solving an additional equation, 
which connects the on-shell scalar coefficients of the k- and t-matrices. In 
the second step care is required, because only five spin-momentum operators 
are independent and not all choices are numerically safe. We show a very 
good agreement among all the considered methods for the selected neutron- 
neutron and neutron-proton observables. In these calculations we use the 
neutron-proton version of a chiral NNLO potential. In the case of the iter- 
ation method we demonstrate also a fast convergence with respect to the 
number of iterations. This iterative method is best suited to be used in the 
3N bound state and in the 3N continuum calculations. In fact it has been 
already successfully used in the 3N bound state calculations [TS], realizing 
the theoretical formulation given in [T3] . 

Last not least, we apply our 3D framework to the problem of proton- 
proton scattering, using the proton-proton version of the chiral NNLO NN 
potential, supplemented with the (screened) Coulomb force. We show con- 
verged results for two different screening functions and find a good agreement 
with two other methods considered in [27] . 
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Fig. 1 Selected observables for the neutron-neutron (left panel) and neutron- 
proton (right panel) system at the projectile laboratory kinetic energy 13 MeV 
as a function of the center-of-mass angle 6 c . m . for the chiral NNLO potential [19] - 
Crosses represent results obtained with the iterative method (see text). Dashed 
(dotted) lines depict the results obtained by solving first the equation for the It- 
matrix and then the on-shell equation for the t-matrix, where on shell the W4 (we) 
2N operator is eliminated. Solid lines are for the predictions based on the mat rix 
inversion method. For the definition of the R, A and D observables see e.g. |17] . 
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Fig. 2 The same as in Fig. Q] for the projectile laboratory kinetic energy 
£=300 MeV. 
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Fig. 3 The Wolfenstein parameters for neutron-proton scattering for proje ctile 
laboratory kinetic energy 300 MeV calculated with the chiral NNLO potential 19, . 
The left panels show the real parts of the amplitudes, whereas the imaginary parts 
are displayed in the right panels. Crosses and lines show results of the same types 
of calculations as in Fig. \T\ 
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Fig. 4 The convergence of results obtained with the iterative method using dif- 
ferent numbers of iterations for neutron-neutron (left panel) and neutron-proton 
(right panel) scattering observables for the projectile laboratory kinetic energy of 
13 MeV and the chiral NNLO potential [T5]. Crosses show the predictions with 6 
iterations. Dash-dotted, double-dotted, dotted, dashed and solid lines display the 
results calculated with 8, 10, 12, 14 and 20 iterations, respectively. The last three 
lines overlap. 




Fig. 5 The same as in Fig. [4] for the projectile laboratory kinetic energy of 
300 MeV. 
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Fig. 6 Selected observables for the proton-proton system (including the Coulomb 
force) at the projectile laboratory kinetic energy 13 MeV as a functi on of th e center- 
of-mass angle 6 c . m . for the chiral proton-proton NNLO potential [19125] . Results 
are obtained by solving the LS equation using the matrix inversion method and 
applying the exponential screening, si(r;n, R), to the Coulomb. Dotted, dashed 
and solid lines show the results with the screening parameter R— 20, 60 and 120 
fm, respectively, and n — 4. 
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Fig. 7 The same as in Fig. [6] with the localized (sine function) screening (|29fl , 
applied to the Coulomb force. Dotted, dashed and solid lines show the results with 
the screening parameter R— 10, 30 and 55 fm, respectively. 
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Fig. 8 Selected observables for the proton-proton system (including the Coulomb 
force) at the projectile laboratory kinetic energy 13 MeV as a functi on of th e center- 
of-mass angle 8 c . m . for the chiral proton-proton NNLO potential [19125] . Results 
with the si(r; n = 4, R = 120 fm) screening (dashed lines) are compared with the 
predictions based on the S2(r;R = 55 fm) screening (solid lines). In the inset the 
two screening functions are compared: the dashed (solid) line shows the Si (S2) 
function. 
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Fig. 9 Selected observables for the proton-proton system (including the Coulomb 
force) at the projectile laboratory kinetic energy 13 MeV as a funct ion of t he center- 
of-mass angle 6 .m. for the chiral proton-proton NNLO potential [19l|25j . Fully 3D 
results with the Si(r; n = 4, R — 120 fm) screening (solid lines) are compared with 
the predictions based on two methods used in Ref. [27]: the method combining a 
3D Coulomb matrix (for the same scre eni ng si) with PWD calculations (dashed 
lines) and the Vincent-Phatak method 29 (dotted lines). 
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Fig. 10 The same as in Fig. [5] but the S2(r;R — 55 fm) screening [3] is used in 
the 3D calculations involving the Coulomb force. 



